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ESTIMATING  THE  SCATTERING  COEFFICIENT  OF  THE  OCEAN  SURFACE 


FOR  HIGH-FREQUENCY  OVER-THE-HORIZON  RADAR 


INTRODUCTION 


The  radar  scattering  coefficient  of  the  ocean  surface.  O’® .  is  a  basic  quantity  in  radar  remote 
sensing  of  the  ocean  surface  for  such  purposes  as  sea  state  estimation,  weather  and  storm  tracking, 
the  estimation  of  the  propagation  properties  of  the  atmosphere  and  the  quantitative  aspects  of  the 
detection  and  monitoring  of  targets  on  and  above  the  ocean  surface.  One  of  the  most  important 
uses  of  (r°  is  in  the  real-time  assessment  of  a  surveillance  radar’s  performance  capabilities.  If  (r® 
is  known,  the  ratio  of  the  backscatter  from  the  ocean  surface  to  noise  can  be  used  to  indicate  the 
radar's  performance  indices  Pq  (probability  of  detection)  and  PpA  (probability  of  false  alarm). 

.\lso,  an  estimation  of  cr®  can  aid  in  the  determination  of  the  cross-section  of  a  target  detected  by 
an  High-Frequency  Over-the-Horizon  (HF  OTH)  radar  system.  The  cross-section  of  the  target  is 
computed  using  O’® ,  the  area  of  the  radar  footprint  on  the  ocean  surface  and  the  ratio  of  the  power 
level  of  the  target  to  that  of  the  sea  clutter  in  the  doppler  spectrum. 

In  HF  OTH  radar  systems  the  propagation  path  losses  in  terms  of  amplitude  and  polarization 
rotation  are  not  fully  known.  Thus,  it  is  very  difficult  to  estimate  o®  directly  from  the  radar 
equation.  The  technique  developed  by  Trizna  [l]  is  based  on  the  doppler  characteristics  of  the  first 
and  second  order  radar  scatter  from  the  wind  disturbed  ocean  surface,  allowing  an  estimation  of 
O’®  from  internal  doppler  spectrum  power  ratios.  It  was  developed  through  extensive  averaging  of 
ground  wave  data  obtained  using  relatively  long  coherent  integration  times  (CIT). 

The  operational  efficiency  of  an  HF  OTH  system  used  for  surveillance  does  not.  however,  allow 
the  option  of  extensive  averaging  in  time.  The  purpose  of  the  present  investigation  is  to  extend 
the  technique  to  OTH  sky  wave  propagation  and  apply  it  to  relatively  short  CIT  data  and  shorter 
averaging  times.  The  basic  technique  is  reviewed  in  the  next  Section.  This  is  followed  by  a  dis¬ 
cussion  of  some  of  the  possible  methods  of  applying  the  procedure  to  two  sets  of  HF  OTH  radar  I 

data  collected  on  13  August  1981.  The  results  of  the  various  methods  are  reviewed  and  it  is  con¬ 
cluded  that  one  specific  method  yields  the  best  results.  A  method  of  operation  is  recommended 
for  actual  data  collection.  • 

k 

I 
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THE  BASIC  TECHNIQUE  i 


The  technique  developed  by  Trizna  is  based  on  data  collected  using  vertically  polarized  anten¬ 
nas  and  ground  wave  propagation.  It  is  directly  applicable  to  OTH  skywave  propagation  in  yield¬ 
ing  a  scattering  coefficient  identical  to  that  which  would  be  measured  using  vertically  polarized 
antennas  and  ground  wave  propagation,  i.e.  or,.® .  However,  a  3  dB  modification  must  be  introduced 
in  determining  the  cross-section  of  a  target. 

Manuicript  approved  January  7,  1986. 
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The  effective  scattering  coefficient  for  an  OTH  radar  is  smaller  than  (Tvv-  During  propagation 
through  the  ionosphere  radio  waves  experience  polarization  rotation  (Faraday  rotation).  The  effect 
is  such  that  the  polarization  of  the  radar  energy  incident  on  the  ocean  surface  undergoes  several 
rotations  (progressing  mainly  in  range)  within  the  resolution  cell  of  a  single  beam.  Since  the 
scattering  coefficient  for  horizontal  polarization  (at  HF  frequencies  and  near  grazing  incidence)  is 
very  small  compared  to  that  for  vertical  polarization,  the  effective  scattering  coefficient  for  a  radar 
resolution  cell  illuminated  via  skywave  propagation  is  3  dB  smaller  than  o\.v. 

The  usual  target  is.  however,  of  small  spatial  extent  relative  to  the  spatial  period  of  the  rotat¬ 
ing  polarization  of  the  incident  radar  energy.  At  any  given  instant  of  time  the  cross-section  of  the 
target  is  not  being  spatially  averaged  in  the  same  manner  as  that  of  the  ocean  surface.  Therefore. 
Trizna's  technique  for  determining  the  cross-section  of  a  target  on  the  ocean  surface  may  be  ap¬ 
plied  directly  to  OTH  radar  data  if  EITHER  3  dB  is  subtracted  from  OR  the  detected  power 
level  of  the  sea  clutter  return  is  increased  by  3  dB  in  comparing  it  with  that  of  the  target.  In 
determining  the  cross-section  of  a  target  above  the  ocean  surface,  multipath  effects  must  also  be 
considered. 

The  pattern  of  rotating  polarization  incident  of  the  ocean  surface  and  targets  within  the  reso¬ 
lution  cell  also  moves  over  the  surface  with  time.  The  polarization  of  the  radar  energy  incident  on 
a  target  at  any  given  instant  of  time  or  during  a  given  QT  is  unknown.  Therefore,  the  cross- 
section  of  a  target  computed  using  the  data  for  a  given  CIT  is  an  average  cross-section  for  the  unk¬ 
nown  range  of  polarization  vector  orientations  incident  on  the  target  during  the  CIT.  Normally, 
time  averaging  over  a  period  of  one  or  significantly  more  cycles  of  polarization  rotation  is  required 
to  obtain  a  reasonable  estimate  of  the  targets  cross-section.  Thus,  such  cross-section  determina¬ 
tions  are  an  average  over  all  polarizations.  By  virtue  of  its  motion,  an  aircraft  would  probably 
pass  through  many  more  cycles  of  polarization  rotation  than  a  ship  on  station. 

It  should  be  noted  that  the  values  given  for  the  various  estimates  for  O’*  in  this  report  ARE 
NOT  EFFECTIVE  SCATTERING  COEFFICIENTS  unless  otherwise  stated.  They  are  instead  esti¬ 
mates  of  OV-v- 

.4n  example  of  an  HF  OTH  radar  doppler  spectrum  obtained  for  a  CIT  of  51.2  seconds  is 
shown  in  Fig.  1.  The  two  first  order  Bragg  peaks  and  the  second  order  return  appear  near  zero  Hz 
doppler  frequency.  The  return  at  +1.1  Hz  is  from  a  ship.  The  estimation  of  O’®  relies  on  the 
measures  rho.  the  power  ratio  of  the  two  first  order  Bragg  peaks,  and  zeta,  the  power  ratio  of  the 
stronger  first  order  Bragg  peak  to  the  second  order  power  level  midway  between  the  two  first  ord¬ 
er  Bragg  peaks.  The  radar  cross-section  of  the  ship  may  then  be  computed  using  <T° .  the  area  of 
the  ocean  surface  illuminated  by  the  radar  beam  and  the  power  ratio  mu  of  the  target  to  the  total 
sea  clutter.  In  this  instance  the  power  level  of  the  total  sea  clutter  return  may  be  approximated 
by  that  of  the  stronger  first  order  Bragg  peak. 

The  radar  scattering  coefficient  of  the  ocean  surface  is  determined  using  the  nomograph  shown 
in  Fig.  2  and  the  equation  [l] 

cr®  =  -  38.9  +  A  +  (p/2)  +  10  log  [l  +  10"'’'*®]  . 

where  <r® ,  A  and  p  are  expressed  in  decibels.  In  Fig.  1  rho  *  20.2  dB  and  zeta  -  26.8  dB.  Using 
the  measure  rho,  the  nomograph  indicates  that  the  acute  angle  theta  between  the  azimuthal  direc¬ 
tion  of  propagation  of  the  radar  energy  and  the  major  axis  of  the  ocean  wave  spectrum  was  38  de¬ 
grees.  The  A  -  0  dB  reference  line  for  a  CIT  of  51.2  seconds  is  0.1  dB  below  that  shown  for  CIT  - 
50  seconds  and  a  »  0.(X)81  (Phillips'  constant).  The  value  for  A  to  be  used  in  the  above  equation 
is  seen  to  be  +2.4  dB.  Therefore,  O’®  -  -26.4  dB(sqm/sqm). 
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Fig.  1.  .4n  example  of  HF  OTH  radar  data  displaying  the  received  power  as  a  function  of 
doppler  frequency. 


Fig.  2.  The  nomograph  used  in 
determining  (7°  (adapted  from 
ref.  [1]). 
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The  effective  scattering  coefficient  of  the  resolution  cell  was  then  <r°  -  3.0  -  -29.4 
dB(sqm/sqm).  Since  its  area  was  90.6  dB(sqm),  the  effective  cross-section  of  the  surface  within 
the  resolution  cell  was  61.2  dB(sqin).  The  cross-section  of  the  target  relative  to  the  effective 
cross-section  of  the  sea  clutter  was  the  measure  mu  (-40.3  dB).  yielding  a  cross-section  for  the 
ship  of  20.9  dBCsqm).  This  value  for  the  cross-section  of  the  ship  is  for  the  unknown  range  of 
polarization  vector  orientations  incident  on  the  ship  during  the  51.2  second  CIT.  Based  on  t^e 
computations  reviewed  in  the  next  Section,  the  radar  scattering  coefficient  of  the  ocean  surface,  the 
effective  cross-section  of  the  surface  and  the  cross-section  of  the  ship  indicated  by  the  data  from 
this  particular  CIT  are  too  large  by  about  5  dB.  This  is  an  example  of  one  of  the  extremes  encoun¬ 
tered  in  the  data. 

For  data  analysis  purposes,  equations  have  been  used  in  place  of  the  nomograph.  The  ocean 
wave  spectral  spreading  function  assumed  by  Trizna  determined  the  upper  curve.  Its  functional 
form  is  given  in  Fig.  2.  The  lower  empirical  curve  has  been  approximated  for  this  study.  The  er¬ 
ror  in  the  computed  value  for  A  (a  =  0.0081)  with  respect  to  the  nomograph  is  less  than  ±  0.4  dB 
for  0  ^  ^  theta  <  45  ^  and  +  0.0.  -  0.1  dB  for  45  “  <  theta  <  90  °  . 

This  basic  technique  for  estimating  O’"  is  dependent  on  the  ocean  wave  spectrum  of  the  water 
wavelengths  being  sensed  by  the  radar  [2].  Oceanographic  knowledge  of  the  ocean  wave  spectrum 
is.  however,  incomplete.  The  areas  in  which  the  incompleteness  of  oceanographic  knowledge 
affects  the  determination  of  (r°  are  the  azimuthal  shape  of  the  spectrum  as  a  function  of  water 
wavelength,  and  hence  radar  wavelength,  and  the  shape  of  the  spectrum  with  fetch  and  with  time, 
as  discussed  below. 

If  a  steady  wind  has  been  blowing  over  a  period  of  several  hours,  and  over  a  long  expanse  (a 
few  hundreds  of  miles),  or  fetch,  of  ocean  surface  toward  the  area  being  measured,  the  azimuthal 
spreading  factor  or  shape  of  the  spectrum  for  water  wavelengths  of  about  10  m  and  shorter  is  ap¬ 
proximately  cosine  squared,  as  assumed  here.  The  quantitative  manner  in  which  the  spreading 
factor  differs  from  cosine  squared  for  longer  water  wavelengths  and  shorter  fetches  and  times  is 
not  completely  known. 

For  areas  of  the  ocean  surface  towards  which  the  wind  has  been  blowing  steadily  for  at  least  3 
to  6  hours,  and  over  a  fetch  of  at  least  150  nautical  miles,  the  estimated  error  in  o'" .  when  sensed 
at  radar  frequencies  greater  than  15  MHz.  is  ±  3  dB.  Water  wavelengths  of  10  m  and  shorter 
respond  fairly  rapidly  to  surface  winds,  and  thus  this  region  of  the  ocean  spectrum  can  be  as¬ 
sumed  to  be  in  equilibrium  for  the  above  conditions.  For  radar  frequencies  of  5  to  15  MHz  the  es¬ 
timated  error  is  ±  6  dB.  and  ±  10  dB  for  frequencies  of  1  to  5  MHz.  These  errors  increase  as  the 
angle  between  the  direction  of  the  wind  and  the  direction  of  propagation  of  the  radar  energy  ap¬ 
proaches  90  °  .  i.e.  as  rho  approaches  0  dB. 

This  method  is  assumed  to  be  valid,  but  with  increased  error,  for  shorter  periods  of  a  steady 
wind  after  near  calm  conditions  and  for  a  short  period  after  the  weather  system  starts  to  decay. 
This  method  is  not.  however,  applicable  to  areas  experiencing  changing  meteorological  conditions. 
Therefore,  the  applicability  of  this  method  of  determining  tr"  can  be  enhanced  through  the  use  of 
surface  meteorology  forecasts  and  hindcasts. 
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It  should  be  noted  that  <r° .  with  or  without  identifying  subscripts,  should  be  read  as  cTyv  with 
or  without  the  same  identifying  subscripts. 
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Fig.  3.  A  falling  raster  display  of 
the  received  power  as  a  function  of 
doppler  frequency  for  successive 
coherent  integration  times  (CIT  = 
51.2  s)  in  the  interval  16:28:30  to 
16:38:02  GMT  on  13  August  1981. 
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The  falling  raster  display  shown  in  Fig.  3  is  a  plot  of  the  received  power  as  a  function  of 
doppler  frequency  for  successive  coherent  integration  times.  All  of  the  returns  are  from  the  same 
range,  as  indicated  by  the  range  in  nautical  miles  along  the  right  edge  of  the  Figure,  just  above  the 
time  (GMT)  for  each  CIT.  The  termini  of  the  lines  connecting  the  center  frequencies  of  the  two 
first  order  Bragg  peaks  indicate  the  relative  powers  of  the  two  peaks  based  on  the  median  power 
levels  of  similar  data.  The  vertical  separation  between  the  lines  is  a  constant  25  dB  -  the  power 
offset  in  the  falling  raster  display.  The  center  points  of  the  short  lines  through  the  second  order 
Bragg  return,  midway  between  the  two  first  order  peaks,  also  indicate  the  relative  power  levels  of 
similar  data.  The  purpose  of  these  lines  is  to  represent  reference  power  levels  against  which  CIT 
to  CIT  variations  in  power  may  be  visualized.  It  is  seen  that  the  power  levels  of  interest  in  deter¬ 
mining  rho  and  zeta  vary  both  in  absolute  and  relative  magnitude. 
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Fig.  4.  A  falling  raster  display  of  the  received  power  as  a  function  of  doppler  frequency  for 
successive  coherent  integration  times  (CIT  «  12.8  s)  in  the  interval  16:38:40  to  16:44:33 
GMT  on  13  August  1981. 
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Fig.  5.  The  average  power,  at  each  dopplet  frequency,  for  the  eleven  traces  shown  in  Fig.  3. 


The  data  shown  in  Fig.  3  were  obtained  with  a  CTT  of  51.2  seconds.  Data  obtained  immediate¬ 
ly  after  that  displayed  in  Fig.  3.  but  with  a  CIT  of  12.8  seconds,  are  similarly  displayed  in  Fig.  4. 
These  data,  as  well  as  other  data  sets,  also  indicate  that  the  power  levels  of  interest  in  determining 
rho  and  zeta  vary  erratically  from  CIT  to  CTT. 

Fig.  5  displays  the  average  power  (the  square  of  the  average  of  the  square  root  of  the  power 
levels)  at  each  doppler  frequency  for  the  eleven  traces  shown  in  Fig.  3.  The  reference  power  levels 
are  the  same  as  those  shown  in  Fig.  3.  The  values  of  rho  and  zeta  measured  from  the  average 
power  levels  are  15.5  and  25.2  dB,  yielding  a  cr®  of  -31.2  dB(sqm/sqm).  .Averages  of  the  first, 
second  and  third  sets  of  nine  traces  from  the  27  shown  in  Fig.  4.  and  the  average  of  all  27.  are 
similarly  displayed  in  Figs.  6.  7.  8  and  9.  The  values  of  rho  and  zeta  measured  from  the  average 
power  levels  yield  values  for  <r°  of  -35.6.  -38.4.  -33.2  and  -35.5  dB(sqm/sqm).  Thus,  the  average 
power  levels  for  these  data  sets  also  exhibit  significant  disparities  in  the  value  for  <r® . 

In  computing  the  cross-section  of  a  moving  target,  either  a  ship  or  an  aircraft,  the  scattering 
coefficient  would  be  determined  for  the  range  bins  which  follow  the  track  of  the  target.  The  re¬ 
turn  from  a  ship  is  apparent  in  Figs.  3  and  4  at  a  doppler  frequency  of  approximately  +1.1  Hz. 
The  ship  was  tracked  through  both  time  periods  bas^  on  its  maximum  power  return.  Of  the 
eleven  coherent  integration  times  in  the  interval  from  1628  to  1638  hours  (Fig.  3).  the  ship  was 
determined  to  be  in  range  bin  28  (968.7  nautical  miles)  seven  times  and  in  range  bin  27  (960.6 
nmi)  a  total  of  four  times,  as  it  slowly  approached  the  radar.  Of  the  27  coherent  integration  times 
in  the  interval  from  1638  to  1645  hours  (Fig.  4).  the  ship  appeared  in  range  bin  27  a  total  of  15 
times,  with  a  small  number  of  extreme  excursions  from  range  bins  24  (936.3  nmi)  to  30  (993.0 
nmi). 
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Fig.  6.  The  average  power,  ai  each  doppler  frequency,  for  the  first  set  of  nine  traces  shown 
in  Fig.  4. 
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Fig.  7.  The  average  power,  at  each  doppler  frequency,  for  the  second  set  of  nine  traces 
shown  in  Fig.  4. 


Fig.  10.  The  variation  in  the 
received  power  levels  (using  the 
peak  width  criteria 

MXPK\VD=1 )  for  successive 
coherent  integration  times,  and 
the  resulting  values  for  <T* 
computed  for  each  CIT. 
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The  data  obtained  from  tracking  the  ship  have  been  processed  in  three  different  ways.  The 
results  are  displayed  in  the  lower  halves  of  Figs.  10.  11  and  12.  In  Fig.  10,  MXPKWD,  the  max¬ 
imum  peak  width  input  for  the  tracking  programs,  is  unity.  Thus,  the  power  levels  plotted  for 
the  recede  and  approach  peaks  and  the  central  peak  are  the  maximum  power  levels  attained  in  a 
single,  discrete  doppler  frequency  bin.  The  recede  peak  is  the  first  order  Bragg  peak  due  to  waves 
propagating  away  from  the  radar  and  the  approach  peak  is  that  due  to  waves  propagating  towards 
the  radar.  A  central  peak  is  a  peak  appearing  within  the  central  20%  (or  5  doppler  bins,  whichev¬ 
er  is  greater)  of  the  doppler  interval  between  the  two  first  order  Bragg  peaks.  For  the  51.2  second 
data,  in  which  the  separation  of  the  two  first  order  Bragg  peaks  ranged  from  45  to  47  doppler  bins 
(0.0195  Hz  per  bin),  the  central  region  is  the  middle  9  doppler  bins  (0.176  Hz).  The  average 
separation  for  the  12.8  second  data  was  nearly  11.6  doppler  bins  (0.0781  Hz  per  bin).  For  such 
data  the  central  20%  would  be  too  small  an  interval  to  search  for  a  second  order  peak.  Thus,  the 
central  region  was  defined  to  be  the  middle  5  doppler  bins  (0.391  Hz).  A  central  peak  cannot  al¬ 
ways  be  found.  One  is  missing  in  the  data  from  1628  to  1638  hours  and  only  eight  were  found  in 
the  data  for  1638  to  1645  hours.  The  mid  doppler  power  level  is  that  in  the  doppler  bin  midway 
between  the  doppler  bins  of  the  two  first  order  Bragg  peaks. 

In  Fig.  11.  MXPKWD  is  three.  Thus,  the  power  levels  were  obtained  by  integrating  the  square 
root  of  the  powers  in  the  doppler  bin  of  maximum  power  (the  mid  doppler  bin  in  the  case  of  mid 
doppler)  and  one  doppler  bin  on  each  side  of  it.  The  integrated  amplitudes  were  then  squared  to 
yield  the  appropriate  power  levels.  In  Fig.  12.  where  MXPKWD  is  19,  the  integration  is  over  the 
doppler  bin  of  maximum  power  return  ±  9  bins  or  the  full  width  of  the  peak,  whichever  is  less. 
The  full  width  of  a  peak  is  taken  to  be  the  doppler  width  between  the  points  at  which  the  powers 
increase  to  form  adjacent  peaks.  The  doppler  interval  resulting  from  integrating  over  the  mid 
doppler  bin  ±  9  bins  would  be  too  broad.  Therefore,  that  data  is  not  included  in  Fig.  12.  The  er¬ 
ratic  fluctuations  in  the  power  levels  evident  in  Figs.  3  and  4  are  also  evident  in  the  lower  halves 
of  Figs.  10,  11  and  12. 
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Fig.  11.  The  variation  in  the 
received  power  levels  (using  the 
peak.  width  criteria 

■MXPKWDmJ)  for  successive 
coherent  integration  times,  and 
the  resulting  values  for  <r° 
computed  f  jr  each  CIT. 


-  -A  A' 

•  •  ‘.e.  .  V, 


''  A  )  A 

tV-V'  ' 


y.. 


—  —  4PM0ACH  PEAK 


Fig.  12.  The  variation  in  the 
received  power  levels  (using  the 
peak  width  criteria 

MXPKWD—19)  for  successive 
coherent  integration  times,  and 
the  resulting  values  for  (r° 
computed  for  each  CTT. 
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These  three  different  ways  of  measuring  the  power  in  a  peak  are  being  investigated  in  an  effort 
to  determine  the  measure  that  will  yield  the  best  practical  results.  The  first  measure,  MXPKWD 
“  1,  is  the  simplest  to  implement.  However,  if  the  power  in  the  mid  doppler  bin  is  used  to  deter¬ 
mine  zeta,  that  power  level  may  be  at  the  peak  of  a  central  peak,  in  the  trough  between  two  near¬ 
central  peaks  or  at  some  power  level  between  these  two.  Thus,  zeta  may  be  too  large,  leading  to  a 
low  estimate  for  O’*' .  The  use  of  MXPKWD  *  3  or  slightly  larger  may  alleviate  this  problem.  .A 
maximum  peak  width  of  19  has  been  used  for  tracking  targets.  A  cursory  examination  of  some  of 
the  available  data  indicated  that  the  full  peak  width  of  most  of  the  targets  was  about  eleven 
doppler  bins.  Nineteen  was  chosen  to  accommodate  some  peak  broadening.  The  full  width  of  a 
central  peak  is  usually  much  less  than  19  doppler  bins.  The  two  first  order  peaks  are.  however, 
much  wider.  This  disparity  may  also  lead  to  a  high  value  for  zeta  and  a  low  estimate  for  <T° . 

An  estimation  of  <T°  requires  a  determination  or  estimation  of  the  measures  rho  and  zeta. 
While  it  can  be  seen  from  Figs.  3.  4.  10,  11  and  12  that  the  value  for  rho  fluctuates  erratically,  the 
measurement  itself  is  straightforward.  However,  the  near  continuum  of  the  second  order  scatter 
midway  between  the  two  first  order  Bragg  peaks  that  is  available  when  using  relative  long 
coherent  integration  times  and  extensive  averaging,  is  not  evident  in  the  data  displayed  here. 
Thus,  what  second  order  power  level  should  be  used  to  determine  zeta? 

Three  different  second  order  power  levels  will  be  investigated.  They  are  the  power  level  of  a 
central  peak,  the  power  level  in  the  doppler  bin  midway  between  the  two  first  order  Bragg  peaks, 
and  a  power  level  based  on  the  measure  rho  and  the  assumption  of  ocean  wave  spectral  saturation 
( A  »  0)  at  the  appropriate  water  wavelengths. 

The  various  ways  in  which  <r°  may  be  estimated  from  these  power  levels  are  indicated  in  Fig. 
13.  as  applied  to  a  sample  of  the  data  extracted  from  Fig.  11.  The  estimates  <t^  and  it^°  are  based 
on  the  power  levels  of  the  central  peaks.  The  estimate  (r^  is  determined  by  first  computing  a  cr„° 
from  the  recede,  approach  and  central  peak  power  levels  for  each  CTT.  A  mean  and  median  are 
then  computed  from  the  <r,®  values.  .Alternatively,  is  determined  by  first  computing  the  mean 
and  median  of  the  dB  power  levels  for  the  recede,  approach  and  central  peaks  in  the  data  set,  and 
then  computing  (mean  and  median).  The  estimates  di®  and  (T^  are  determined  in  a  similar 
manner,  using  the  mid  doppler  power  levels  in  place  of  the  central  peak  power  levels.  The  esti¬ 
mates  O’e”  and  crf°  are  based  on  the  recede  and  approach  peak  power  levels  and  the  assumption  of 
ocean  wave  spectral  saturation.  The  estimate  <7^°  is  determined  by  first  computing  a  <r„°  using  the 
measure  p„  and  A  -  0  for  each  CIT.  A  mean  and  median  <r^°  are  then  computed  from  the  <t„° 
values.  Alternatively,  tr/  is  determined  by  first  computing  the  mean  and  median  of  the  dB  power 
levels  for  the  recede  and  approach  peaks  in  the  data  set  and  then  computing  p  and  <7^°  (mean  and 
median)  assuming  A  »  0.  These  six  methods  of  estimating  <7°  will  be  discussed  below. 

The  values  of  or,®  computed  from  the  data  available  in  each  CIT  are  displayed  in  the  upper 
halves  of  Figs.  10,  1 1  and  12.  The  values  based  on  the  central  peaks  are  indicated  by  the  black¬ 
ened  circles.  The  mean,  standard  deviation  and  median  of  the  values  for  each  time  period, 
identified  by  the  coherent  integration  time  and  MXPKWD,  are  tabulated  in  TABLE  1.  There  are 
too  few  data  points  in  most  of  these  data  sets  to  allow  a  meaningful  computation  of  the  standard 
deviation.  These  values  will  only  be  utilized  for  a  rough  comparison  of  the  spread  in  values  of  (r^® 
and  Og®  with  those  of  erg®. 

The  experimental  values  for  o'®  may  be  compared  against  a  theoretical  value  for  <7°  based  on  a 
cosine  squared  spreading  function  for  the  ocean  spectrum  and  a  total  spectral  energy  content  in 
agreement  with  the  work  of  Barrick  [3],  as  in  ref.  [l].  It  will  be  assumed  that  the  ocean  spectrum 
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Fig.  13.  sample  of  the  data  extracted  from  Fig.  11  and  some  of  the  various  vays  in  vhich 
an  estimate  of  O’"  may  be  determined  using  the  recede  peak  (R),  approach  peak  (A),  central 
peak  (C)  or  mid  doppler  (M)  power  levels  for  each  CIT,  or  R^.  A^.  Cm  or  Mm.  tbe  mean 
(or  median)  of  the  dB  power  levels  for  the  data  set. 


at  the  water  wavelengths  scattering  the  19.071  MHz  radar  energy,  i.e.  7.9  meters  for  the  two  first 
order  Bragg  peaks  and  5.6  meters  for  the  second  order  return  at  mid  doppler,  was  saturated.  The 
mean  of  the  dB  values  of  rho  for  all  111  pairs  of  first  order  Bragg  peaks  plotted  in  Figs.  10.  11 
and  12  is  14.3  dB.  and  the  median  is  14  dB.  This  implies  an  angle  theta  of  58  ^  and  a  theoretical 
value  for  tr®  of  -31.6  dB(sqm/sqm).  This  value  is  less  than  ten  of  the  twelve  values  for  <r^®. 
lending  credence  to  the  assumption  of  spectral  saturation  at  7.9  and  5.6  meter  water  wavelengths. 

While  the  values  of  <r^®  for  a  CIT  of  51.2  seconds  is  in  good  agreement  with  the  theoretical 
value,  the  values  for  12.8  seconds  are  mostly  several  dB  too  high. 


The  values  of  tr„°  based  on  the  mid  doppler  power  levels  are  joined  by  dotted  lines  in  the  upper 
halves  of  Figs.  10  and  11.  The  mean,  standard  deviation  and  median  values  of  (r,®  are  also  tabu¬ 
lated  in  TABLE  1.  It  is  seen  that  these  values  are  consistently  lower  than  those  based  on  the  cen¬ 
tral  peaks.  The  standard  deviations  are.  however,  comparable. 


Values  for  <7®  may  also  be  computed  from  the  mean  and  median  of  the  dB  power  levels  for 
each  of  the  six  sets  of  data  displayed  in  the  lower  halves  of  Figs.  10.  11  and  12.  The  values  com¬ 
puted  from  the  recede,  approach  and  central  peak  power  levels  are  tabulated  in  TABLE  1  under 
the  Cc  heading.  It  is  seen  that  these  values  are  comparable  to  those  for  for  a  CIT  of  51.2 
seconds  and  that  the  values  for  a  CIT  of  12.8  seconds  are  now  comparable  to  those  for  51.2 
seconds  and  the  theoretical  value  of  -31.6  dB(sq m/sqm). 
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The  values  for  0"“  compuied  from  the  mean  and  median  of  ihe  dB  power  levels  of  ihe  recede 
and  approach  peaks  and  mid  doppler  are  tabulated  in  TABLE  1  under  the  heading  Like  those 
for  (T^'\  these  values  are  consistently  lower  than  those  based  on  the  power  levels  of  the  central 
peaks. 


The  theoretical  value  for  cr°  of  -31.6  dB(sqm/sqm)  was  computed  assuming  that  the  ocean 
spectrum  at  vater  wavelengths  of  7.9  and  5.6  meters  was  saturated.  This  is  a  reasonable  assump¬ 
tion  at  these  wavelengths  and  the  experimental  values  compuied  thus  far  indicate  that  the  as¬ 
sumption  was  valid  for  these  sets  of  data.  If  the  same  assumption  is  made  in  computing  O’,®  for 
for  each  CIT.  only  rho  has  to  be  determined.  The  values  of  rho  for  the  data  displayed  in  Figs.  10. 
11  and  12  are  displayed  in  Fig.  14.  The  corresponding  values  of  <T^°  (A  =  0)  are  shown  along  the 
right  edge  of  the  Figure  for  the  theoretically  valid  interval  0  dB  5  p  «  24  dB.  .As  indicated  above, 
the  mean  of  all  111  values  for  rho  is  14.3  dB  and  the  median  is  14  dB.  These  two  measures  yield 
an  average  value  for  <7°  of  -31.6  dB(sqm/sqm).  The  mean,  standard  deviation  and  median  values 
for  O’®  computed  directly  from  the  values  of  rho  (A  =  0)  for  each  individual  data  set  are  tabulated 
under  the  heading  o’^®  in  TABLE  1.  The  mean  and  median  values  are  of  course  very  close  to  the 
average  value  of  -31.6  dB(sq m/sqm). 

The  standard  deviation  of  the  values  are  mainly  one-half  to  one-third  of  those  tabulated  for 
o’^®.  This  indicates  that  in  computing  o’®  for  each  CIT  using  rho  and  zeia.  one-third  to  one-half  of 
the  spread  in  the  value  of  <T^  is  due  to  the  variability  in  the  measurement  of  rho  and  one- half  to 
two-thirds  of  the  spread  in  due  to  the  variability  in  the  measurement  of  zeta.  This  also  holds  for 
a  comparison  of  the  standard  deviations  for  O’g®  and  cr,®. 


The  assumption  of  ocean  spectrum  saturation  at  7.9  and  5.6  meter  water  wavelengths  may 
also  be  applied  to  the  mean  and  median  of  the  dB  power  levels  for  each  data  set.  The  results  are 
tabulated  under  the  heading  <7^°  in  TABLE  1.  These  values  are  in  good  agreement  with  the 
theoretical  value  for  <7°  of  -31.6  dB(sqm/sqm). 


The  assumption  of  spectral  saturation  yields  a  "maximum"  value  for  O’®.  If  O’®  is  computed 
in  the  manner  of  either  o’/  or  a/,  the  actual  value  of  O’®  would  be  lower  by  the  amount  that  the 
ocean  spectrum  is  below  saturation  at  the  water  wavelengths  responsible  for  the  first  and  second 
order  scattering  of  the  radar  energy. 


The  estimates  O’/,  o’/  and  O’p®  have  been  computed  from  the  mean  and  median  of  the  dB  power 
levels  for  ensembles  of  coherent  integration  times.  These  mean  and  median  levels  may  be  ob¬ 
tained  from  probability  plots,  such  as  those  shown  in  Fig.  15.  These  computer  generated  plots  are 
labeled  in  units  of  the  standard  deviation  sigma.  The  toul  plot  error  is  less  than  ±  1  division  ex¬ 
cept  in  the  six  columns  along  each  side  of  the  plot,  where  the  error  increases  rapidly. 


The  recede  and  approach  power  levels  are  identified  by  "R"  and  ".A",  central  peaks  by  "Z"  and 
mid  doppler  power  levels  by  "M".  The  power  return  from  the  ship  is  identified  by  *T*  and  the 
"+"  and  symbols  represent  the  noise  power  levels  in  the  extreme  positive  and  negative  doppler 
bins. 
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TABLE  1 


ir/iay 
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mcdiftti  mean  mtdian  mean  mtdian  mean  medtan  mean  median 


■U.'y  38  30 

31. b  3.0  30.3 

32.  f  43  32 


-395  6.4 

35.9  .  2.3 


28.3  5.9  -25  5  .  -35.7  6-5  -35  30.4 

28.4  6.3  -27  -34.3  5-9  -32-5  31.0 

31.1  77  -29  -  —  -  32.1 


31-6  1.9  -32.2  31  7  -32.2 

-31.7  2.0  -31.7  ;  -31.8  -32.4 

31  8  1.3  -31.7  31.8  ;  -31. 7  ^ 

- - - ^ ^ - 

-31.0  ;  2.6  -31.0  j  -31.0  -30.8  ^ 

-31.4  2.5  -31.7  I  -31.6  :  -31.7 

-31.7  ■  2.2  -31.7  j  -31.8  ;  -31.7 


Computed  using  the  recede  peak,  approach  peak  and  central  peak  power  level  for  each  CIT. 

Computed  using  the  recede  peak,  approach  peak  and  mid  doppler  power  level  for  each  CIT. 

Computed  using  the  mean  ( median)  of  the  dB  po^i^er  levels  for  the  recede  peaks,  approach  peaks  and  central  peaks  m 
each  data  set. 

Computed  using  the  mean  (median)  of  the  dB  pow’er  levels  for  the  recede  peaks,  approadt  peaks  and  mid  doppler  fre> 
quency  in  each  data  set. 

Computed  using  the  recede  and  approach  peak  power  levels  for  each  CIT  (i^O). 

Computed  using  the  mean  (median)  of  thedB  power  levels  for  the  recede  and  approach  peaks  in  each  data  set  (^0). 
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Fig.  14.  The  values  of  rho  for  all  of  the  data  displayed  in  Figs.  10,  11  and  12  (where 
MKEKWO  «  1,  3  and  19)  and  the  corresponding  values  for  ,  assuming  spectral  satura¬ 
tion  at  water  wavelengths  of  7.9  and  S.6  meters  ( i.e.  A  >  0). 
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The  data  for  each  power  level  are  seen  to  approximate  straight  lines,  indicating  a  Gaussian  dis¬ 
tribution.  The  median  may  be  taken  as  the  power  level  at  sigma  =  0  or  the  average  of  the  two  or 
three  power  levels  nearest  sigma  =  0.  A  straight  line  may  also  be  fitted  to  the  data  between  sigma 
equals  -1  and  +1.  -2  and  +2  or  some  other  limits.  The  mean  power  level  would  then  be  taken 
where  the  straight  line  intersects  sigma  =  0. 
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Fig.  15.  Probability  plots,  in  units  of  the  standard  deviation  sigma,  of  the  power  levels 
(MXPKVD-3)  for  the  recede  (R),  approach  (.\),  central  peak  (Z),  mid  doppler  (M)  and  ship 
(T)  returns  and  the  noise  power  levels  (MXPKWIVl)  in  the  extreme  positive  (•*•)  and  nega¬ 
tive  doppler  bins. 


REVIEW  OF  RESULTS 


The  various  methods  for  estimating  cr"  that  have  been  discussed  in  the  previous  Section  have 
been  applied  to  a  very  small  data  set.  The  major  conclusions  that  may  be  drawn,  however,  are 
consistent  with  the  methodology  and  probably  applicable  in  general. 

The  most  self-consistent  values  for  cr“  have  been  obtained  when  only  the  two  first  order 
Bragg  peaks  are  considered  (  <T^  and  CTp’  ).  These  peaks  represent  the  two  strongest  returns  in  the 
radar  data.  This  fact,  together  with  the  reasonable  assumption  of  saturation  of  the  ocean  wave 
spectrum  at  the  wavelengths  responsible  for  the  radar  return,  lend  credence  to  either  of  these  two 
methods  for  estimating  a-° .  The  spread  in  the  values  for  (mean)  and  crp“  (mean)  is  0.8  dB. 
while  the  spread  is  1.2  dB  for  (median)  and  1.6  dB  for  cTp"  (median).  The  smaller  spread  in 
the  mean  values  with  respect  to  the  median  values  is  also  apparent  in  the  other  methods  for  es¬ 
timating  cr" . 

The  assumption  of  saturation  of  the  ocean  wave  spectrum,  although  reasonable  in  some  in¬ 
stances.  is  not  generally  valid.  Therefore,  the  method  yielding  may  yield  more  reliable  results 
when  used  on  a  continuous  basis  or  as  a  spot  check  on  spectral  saturation. 

The  method  yielding  makes  use  of  the  two  strong  first  order  Bragg  peaks  and  the  central 
peak,  when  it  is  evident.  This  method  yields  the  second  most  consistent  values  for  cr®.  values 
that  are.  in  general,  less  than  1  dB  larger  than  those  for  (r^°  and  cr/  .  The  spreads  in  the  values  for 
<r^°.  1.7  dB  for  the  mean  and  2  dB  for  the  median,  are  greater  than  those  found  for  and  CTp®. 

The  main  problem  with  this  method  is  that  a  central  peak  is  not  always  evident.  One  central 
peak  (9%)  was  missing  in  the  51.2  second  CIT  dau  and  19  (70%)  were  not  evident  in  the  12.8 
second  CIT  data.  Thus,  a  combination  of  the  method  yielding  with  that  for  either  O’/  or  tTp® 
may  produce  the  best  results,  especially  for  short  coherent  integration  times. 

None  of  the  three  methods  reviewed  above  exhibit  a  consistent  variation  with  MXPKWD,  the 
number  of  doppler  bins  over  which  the  energy  in  the  peak  is  integrated. 

The  method  yielding  <r^  is  similar  to  that  for  <r/  except  that  the  power  of  the  central  peak  is 
used  in  the  determination  of  <r„®  for  each  CIT.  It  is  also  similar  to  that  for  tr/  except  that  (Tn®  is 
determined  for  each  CIT  and  then  the  mean  and  median  values  for  O’,®  are  computed,  instead  of 
computing  <r"  from  the  mean  and  median  of  the  dB  power  levels.  The  spreads  in  the  values  for 

are  4.3  dB  for  the  mean  and  6.5  dB  for  the  median.  These  overall  spreads,  and  those  within 
each  subset,  are  greater  than  those  encountered  above.  The  greater  spreads  are  due  to  the  influence 
of  the  individual  central  peak  power  levels.  While  the  values  for  a  CIT  of  51.2  seconds  are  com¬ 
parable  to  those  for  <T^.  {Tp®  and  (Tp®.  five  of  the  six  values  for  a  CIT  of  12.8  seconds  are  several  dB 
larger. 

The  values  for  a,®  and  are  consistently  smaller  than  the  others,  by  as  much  as  6.8  dB. 
This  is  due  to  the  use  of  mid  doppler  power  levels  rather  than  central  peak  power  levels.  The  mid 
doppler  power  level  may  be  at  the  peak  of  a  central  peak,  in  the  trough  between  two  adjacent 
peaks  or  at  some  power  level  between  these  two.  Thus,  zeia  is  too  large  and  cr®  is  too  small. 


CONCLUSIONS 


It  may  be  concluded  from  the  preceding  Section  that  the  best  estimate  of  <T°  is  obtained  using 
the  mean  values  of  the  two  first  order  Bragg  peaks  and  a  central  peak.  Thus,  the  coherent  integra¬ 
tion  time  must  be  long  enough  to  permit  the  detection  and  measurement  of  the  power  in  a  central 
peak. 

The  criteria  for  the  presence  of  a  peak  in  this  work  is  that  the  power  level  two  doppler  bins  on 
one  side  of  the  assumed  peak  must  be  less  than  or  equal  to  the  power  level  of  the  peak  and  that 
the  power  level  two  doppler  bins  on  the  other  side  of  the  assumed  peak  must  be  less  than  that  of 
the  peak.  This  formal  procedure  for  the  establishment  of  a  central  peak  may  be  relaxed  in  estab¬ 
lishing  the  two  first  order  Bragg  peaks.  The  large  amplitude  of  these  two  peaks  relative  to  all  oth¬ 
ers  and  their  approximate  location  and  separation  in  doppler  frequency  indicate  that  these  two 
peaks  may  be  established  by  searching  for  the  maximum  power  levels  at  the  appropriate  doppler 
frequencies. 

The  results  displayed  for  (r^,  Cg'’  and  tr/  in  Table  1  do  not  exhibit  any  consistent  variation 
with  the  degree  of  integration  of  the  energy  in  the  peaks.  Therefore  the  data  processing  may  be 
simplified  by  considering  only  the  peak  power  levels  (MXPKWD  ■  1). 


OPERATIONS  DESIGN  CONSIDERATIONS 


The  best  estimates  of  (r°  are  obtained  for  coherent  integration  times  of  about  30  seconds  and 
longer.  This  is  due  to  the  fact  that  the  variance  of  the  measured  power  levels  decrease  with  time 
and  that  a  central  power  level  becomes  easier  to  establish.  This  conflicts,  however,  with  the  short 
coherent  intergration  time,  or  dwell  time,  required  when  a  large  area  of  the  ocean  surface  is  to  be 
monitored  in  a  timely  manner.  Such  dwell  times  would  be  on  the  order  of  one  second. 

An  integral  part  of  a  short  dwell  time  surveillance  radar  is  an  oblique  sounder  for  frequency 
management  and  performance  capability  assessment.  The  sounder  is  used  to  obtain  ionospheric 
and  other  propagation  information  that  is  not  available  from  short  dwell  time  target  search  data. 
Since  the  revisit  time  constraints  of  the  oblique  sounder  are  not  critical,  it  is  recommended  that  it 
also  be  used  in  a  radar  mode,  with  a  dwell  time  of  30  seconds  or  longer,  to  produce  a  map  of  <r“ 
over  the  area  under  surveillance.  Such  a  map  could  be  used  directly  with  the  short  dwell  time 
search  data  to  establish  the  cross-section  of  detected  targets. 


REFERENCES 

1.  D.B.  Trizna.  "Estimation  of  the  Sea  Surface  Radar  Cross  Section  at  HF  from  Second-Order 
Doppler  Spectrum  Characteristics."  NRL  Report  8579.  May  1982.  AD  A115230 

2.  D.B.  Trizna.  Naval  Research  Laboratory.  Washington.  D.C..  private  communication.  Sept.  1985. 

3.  D.E.  Barrick.  "First-Order  Theory  and  Analysis  of  MF/HF/VHF  Scatter  from  the  Sea,"  IEEE 
Trans.  AP,  AP-20,  pp.  2-10, 1972. 


-  '  *  '  O  O  ■t'V*.*  S*  "AO  S'  ^  *  *'  •.*'  % 


18 


